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In order to describe non-equilibrium media when the manifestations of intrin¬ 
sic structure can not be ignored, we use hydrodynamic mathematical models. 
Information about relaxing processes and interactions between structural ele¬ 
ments is incorporated in the dynamical equations of state (DES) which, unlike 
the local classic relations, now become nonlocal one in time and space. Using 
the symmetry reduction scheme, we obtain the autonomous dynamical systems 
describing the invariant wave solutions. By qualitative analysis methods, we 
show that the dynamical systems possess periodic, multiperiodic, and chaotic 
solutions obeying the Feigenbaum scenario. We study the quasiperiodic regimes 
and their bifurcations. We also reveal the existence of homoclinic trajectories 
of Shilnikov type and investigate the changes of homoclinic structures when 
the bifurcation parameters vary. Hidden attractors, hysteretic phenomena are 
discovered as well. As a result, depending on the chosen model for media, we 
classify the wave solutions and their bifurcations and show that spatio-temporal 
nonlocal models are promising for the describing of complicated wave regimes 
in structured media. 


I. INTRODUCTION 


Open thermodynamic systems attract attention of scientists by their synergetic proper¬ 
ties, their ability to produce localized nontrivial structures and order. Description of such 
phenomena demands the creation of new and the rehnement of already known mathematical 
models. 

According to Refs, lll-y, using the methods of non-equilibrium thermodynamics and the 
internal variables concept^, the nonlinear temporally and spatially nonlocal mathematical 
models for non-equilibrium processes in media with structure have been constructed. In this 
report we present the results of investigations of wave processes in such media. To do this, 
we use the following hydrodynamic type system 
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ctoi cso are the isothermal and adiabatic frozen velocities of sonnd; 700 is the frozen poly¬ 
tropic index. 
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which leads system ([T]) to the dimensionless form 
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We would like to emphasize that system ([3]) can be regarded as an hierarchical set of sub¬ 
models which are complicated by means of taking new effects into account. We thus are 
going to study the chain of nested models and classify their wave solutions using the methods 
of qualitative and numerical analysis. 

The remainder of the report is organized as follows. In SecHIl we begin our studies from 
the simplified version of system ([3]) keeping the terms with the first temporal derivatives, 
then attaching the terms with the second temporal or spatial derivatives. The form of 
wave solutions and the description of techniques for their exploration are presented in de¬ 
tail. Sec lIIIl is devoted to the spatially nonlocal model which is used for investigating of 
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the Shilnikov homoclinic structures whose existence and bifurcations are extremely impor¬ 
tant during chaotic regimes formation. The model incorporated both temporal and spatial 
nonlocalities are presented in Sec lIVI Generalizations of the previous models by means of 
introducing the third temporal derivatives and incorporating of physical nonlinearity are 
given in SecJV] and Sec lVIl respectively. For all models we derive invariant wave solutions 
and carry out the qualitative analysis of corresponding factor-systems. 


II. WAVE SOLUTIONS OF THE MODELS WITH DES INCORPORATING 
THE SECOND TEMPORAL OR SPATIAL DERIVATIVES 


To begin with, let us consider the simplest moc 
5 = h = (T = 0, u = 1 . As has been shown in Refs. \4 

p + pu^ = 0, pu + p^ = 7 p, t{p - xp) = Kp-p, (4) 

due to its symmetry properties^, admits the ansatz 

u = U{oj) + D, p = pq exp(^f -|- 5'(a;)), p = pZ{oj), uj = x — Dt, (5) 


el wit 
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where D is the constant velocity of wave front, ^ determines a slope of the inhomogeneity of 
the steady solution ([5]). According to Ref. fls„>utio„s0a.edescnbedWthep.a„esyste. 
of ODE which possesses limit cycles and homoclinic trajectories. 

If we incorporate the second temporal derivatives in the last equation of system ([3]), then 
the previous DES is generalized to the following one: 
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This model takes into account the dynamics of internal relaxation processes in more detail. 
As has been shown in Ref. Q, wave solutions ([5]) are described by the system of ODE with 
three dimensional phase space. This system possesses the limit cycles undergoing the period 
doubling cascade, and chaotic attractors. 

Consider now the model with relaxation and spatial nonlocality 
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Solutions ([5]) satisfy the following dynamical system 
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which is the only one lying in the physical parameter range. 

We start from analyzing the linearized in the hxed point ([9]) system ([8]) with the matrix 
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The well-known Andronov-Hopf bifurcation theorem^ tells us that periodic solution creation 
can take place if the spectrum of matrix M looks as (—a; ±f2i). This is so if the following 
relations hold: 

a = ^ + 0 0, (10) 

( 11 ) 
( 12 ) 


= AD-B [Zo - D^) + > 0, 

^ {AD - ZqB) > 0. 


The hrst two take on the form of inequalities imposing some restrictions on the parameters. 
The third one determines the neutral stability curve (NSC) in space (H^; n) provided that 
the remaining parameters are hxed. For a = 0.76, ^ = 1.8, r = 0.1, y = 50, it looks like a 
parabola with branches directed from left to right, see Figj2|i. Crossing the NSC from right 
to left, we observe the limit cycle appearance. Development of limit cycle at decreasing D^ 
it is convenient to study by means of the Poincare section techniques'll. 
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Let us choose the plane hh = 0 as an intersecting one and hnd coordinates of intersection 
points of phase curves which cross-sect the intersecting plane only in one direction. Plotting 
coordinate Z of the cross-section point along the vertical axis, and the value of the bifurcation 
parameter along the horizontal one, we will obtain the typical bifurcation diagrams 
(Fig{T]). From the analysis of diagram FigH^ we can see that while parameter decreases 
the development of the limit cycle coincides with the Freihenbaum scenario, followed by the 
creation of a chaotic attractor. Moreover, in the vicinity of the main limit cycle there are 
the hidden attractors (depicted in Fig{T^ by the symbols I and II). These attractors can be 
visualized by the integrating of system ([H]) with special initial data only. 

In Fig{T]D we see the torus development at decreasing D^. According to the diagram, we 
can distinguish tori with densely wound trajectories and striped tori. 

Doing in the same way, we get the two-parameter bifurcation diagram (Fig|2]) which 
tells us that system (IHD possesses the periodic, multiperiodic, quasiperiodic, and chaotic 
trajectories. 

Such a complicated structure of the phase space of the system can be coursed by homo¬ 
clinic trajectory existence. 




FIG. 1. Bifurcation diagrams of system ([8]) in plane (D^, Z), obtained for x = V — 50, C = 1-8, t = 
0.1, a = 0.76 and (a) k = 14, (b) k = 1. 
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FIG. 2. Left: bifurcation diagram of system Q in parametric space 1 - stable focus; 2 - 

IT-cycle; 3 - torus; 4 - multiperiodic attractor; 5 - chaotic attractor; 6 - loss of stability. Right: 
enlargement of part of the left figure: 6 ~ ST-cycle. 


III. HOMOCLINIC LOOPS OF SHILNIKOV TYPE AND THEIR 
BIFURCATIONS 


First it worth noting that existence of honioclinic trajectories, i.e. loops consisting of 
the separatrix orbits of hyperbolic hxed point, plays a crucial role^^i^ in the formation of 
localized regimes in the phase space of dynamical system. 

For the present, the question on the existence of homoclinic trajectory of Shilnikov type^'^ 
in system ((HI) has been treated numerically. 

We investigate a set of points of parameter space (D^, k) for which the trajectories 
moving out of the origin along the one-dimensional unstable invariant manifold LF“ return 
to the origin along the two-dimensional stable invariant manifold LF^. In practice, for the 
given values of par ameters k, D^, we numerically dehne a distance (the counterpart of split 


function in Ref. 
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, p.l98) between the origin and point (A^(a;), Y'"(a;), W^{u)) of the phase 


trajectory F (•; k, D"^): 

f {k, u) = v/iVHi" + [y''(a,)|^ + [W^0J)]\ 
starting from the hxed Cauchy data (0,0,0.001). Next we determine 


$(/?, D^) = min{/'"} 
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FIG. 3. a) Projection of the homoclinic solution of system ([8]) onto the {X, W) plane, b) A portrait 
of subset of parameter space {D‘^, k), corresponding to different intervals of function k) 

values and following Cauchy data; A(0) = y(0) = 0, IF(0) = 0.001: > 1.2 for white colour; 

0-6 < /min < 1-2 for light grey; 0.3 < /foj, < 0-6 for grey; 0.01 < /foj, < 0.3 for deep grey; 
/min ^ 0.01 for black. 


for the part of the trajectory which lies beyond the point at which the distance gains its 
hrst local maximnm, providing that it still lies inside the ball centered at the origin and 
having a hxed (snfficiently large) radins (for this case /^(cu) < 5). The resnlts are presented 
in Figsl3)D. The first is of the most rough scale among this series. Here, white color marks 
the values of parameters k, for which $ > 1.2, light grey corresponds to the cases when 
0.9 < < 1.2 and so on (further explanations are given in the subsequent captions). The 

black coloured patches correspond to the case when <h < 0.01. In Ref. the structure of 
the set of points from Fig. [3b has been studied in more detail. 


IV. MODELS WITH DES TAKING SPATIAL AND TEMPORAL 
NONLOCALITIES INTO ACCOUNT 


Combining the model (E]) and (I7|), we obtain the following spatio-temporal nonlocal model 







































This model has been studied in Refs. [iSl and [la, when the parameters h and a are regarded 
as a small one, i.e., Eqs. (|6]) and ([7]) are perturbed by the terms with high derivatives. It 
turned out that the wave localized regimes are saved under perturbations and undergo some 
smooth changes. 


V. MODELS INVOLVING DES WITH THE THIRD TEMPORAL 
DERIVATIVES 


If we need to describe the relaxing processes in more detail, then we can incorporate the 
terms with the third temporal derivatives in DES fll4p . In this case DES has the form- 
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Solutions ([5]) satisfy the following dynamical system 
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The hxed point of this system has the coordinates 
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where S 3 — — 03 — ag, S 2 — 03 U 9 — a^ag, B\ — 05 ( 03 CI 8 — ciiOg — CI 4 O 7 ), Bq — CI 1 O 5 (0305 — 02 ®?) 
are the coefficients of characteristic polynomial for the matrix M. 

If we hx the parameters x = ?7 = 30, ^ = —1.9, h = 1, r = 1, 6 = 1, a = 2.7, then in the 
plane (S^, k) Eq. (IT^ dehnes the NSC. Crossing this curve in the point A(2.2852; 3.7), one 
can observe the appearance of the limit cycle at D^ > 2.2852. 

In the Poincare diagram depicted at increasing D^ (Fig 01) we can identify the moments 
of several period doubling bifurcations leading to the chaotic attractor creation. But the 
chaotic attractor existing at a short interval of parameter D^ is destroyed. Instead of it 
in the phase space of system flTOl) the complicated periodic trajectory resembling to a loop 
(Fig. [5ti) appears. 

Consider also the development of oscillating regimes whose basins of attraction are sepa¬ 
rated from the basin of attraction of the main limit cycle. Integrating dynamical system flTbl) 
from initial conditions (0; 0; 0; 0.01) at D^ = 2.722, we see that the phase space of the system, 
in addition to the main limit cycle, contains the complicated trajectory (Fig. |5la) which can 
be regarded as a hidden attractor. From the analysis of Poincare diagram (Fig. [ 6 ^) it follows 
that the system weakly responds to the growing of the parameter D^ until D^ = 2.7445. 
When D^ > 2.7445, the system jumps to another type of oscillations followed by chaotic 
regime creation. 

If we plot the Poincare diagram at decreasing D^ (Fig. Eb) starting from the chaotic 
attractor, then we observe the periodic trajectory (FigJSb) that differs from the initial regime 
(FigJSb). Note that the periodic trajectory from Fig. [5b can be revealed directly by the 
integration from the initial conditions ( 0 ; 0 ; 0 ; 0 . 1 ). 
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a) b) 

FIG. 4. a) Neutral stability curve in the plane b) The bifurcation Poincare diagram at 

increasing 




- 0.80 

a) b) 

FIG. 5. Phase portraits of separated trajectories derived at = 2.722, k = 3.7, 5 = 1 and 
different initial conditions. 


VI. DES WITH PHYSICAL NONLINEARITY AND SECOND 
DERIVATIVES 


Till now we dealt with the models without physical nonlinearity. Generalizing the previ¬ 
ous models in this direction, we obtain the following modeP^ 
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a) b) 

FIG. 6. The bifurcation Poincare diagram of development of separated regime at increasing 
(a) and decreasing D^. Here 6 = 1. 

+hxp~‘^[-p (1 + a) + (1 - 0.5a(n - 1)) + Th~^ p (1 + a)] + = p+ (20) 

+Tp - hp - a + p^p^p~^^ , a = (5np”+^ 


Properties of solutions to system fl20ll can be found out using the symmetry of the system 
with respect to the Galilei group>^. One can be persuaded by direct verification that system 
fl20|) allows operator 

- 1 9 d d 

dt ^ ^dx ^ du' 

Let us construct an anzatz with its invariants 


p = R{u), p = P{u), u = 2^t + U{u), u = x — ^R, (21) 

where parameter ^ is proportional to acceleration of the wave front. Substitution fl2Tl) into 
the system yields the following quadrature 

UR = C = const 
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and the dynamical system 
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where (■)' = — (•), L = 1 + a, M = 1 — an, N = 1 — 0.5a(n — 1), a = 6nR^~^^. 
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A=(c^-xL) (C^hRf + a). 

The NSC for system (122|1 has the following form 

G (^, CT, n, h, T, K, y) = 05 (03 T 0204) T 0404 = 0 . 

Let us make the values of parameters hxed as follows: 

7 = 1, x = 10, (7 =-2.8, CT = 0.2, r = 1.1, h = 3.2, 5 = 1.4, n = m = 3.2. 


(26) 
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Condition allows us to find numerically the value of = 0.157 corresponding to 
birth of the limit cycle. 

Let us consider in more detail the influence on the revealed regimes of parameters n and 
5 changes, which determine nonlinearity of the medium in the dynamic equation of state. 
Let us make the value of parameter ^ = 0.35 hxed, in case of which there is a limit cycle 
with period 2T in the space of the system. 

The diagram reveals some peculiarities of system’s 022 p behavior. In particular, we would 
like to pay attention to the presence of a ’’special” point in the parameter plane, surrounded 
by four different types of solutions. One can also see the ’’windows” of periodicity (area 6) 
among the chaotic area. To hnd out the structure of phase space in more detail near area 6 
of FigJTK, let us plot a one-parametric Poincare diagram (FigJTb) for 5 = 1.4 and a decrease 
of parameter n. 

In case of n close to 2, abrupt reconstruction of the chaotic attractor structure can be 
observed, which is probably caused by the interaction of two (or more) co-existing attractors 
of a dynamic system. In case of n ~ 1.4 the chaotic trajectory is localized in a more narrow 
area of phase space of system ( 12 ^ . stipulating the appearance of a specihc window of 
periodicity with a decrease of n. Analysis of a two-parametric bifurcation diagram for the 
value of parameter k = 2 (FigJT^) shows that the area of existence of the chaotic attractor 
increases and the windows of regular intervals in case of the increasing k shift towards higher 
values of the nonlinearity parameter n. 

A crucially different set of bifurcations is observed in case of a change of parameter a. 

Let us £x the values of parameters 7 = 1, y = 10, (7 = —2.8, r = 1.1, k = 0.9, h = 3.2, 
^ = 0.35, n = m = 3.2 and 6 = 0.4. Integrating system fl2^ with initial data (0, 0, 0.01) 
and cr = 5 within phase space near the equilibrium point, in addition to the limit cycle, 
other periodic trajectory has been found with a separated pool of attraction (development 
of this regime with increasing of a is presented in FigJTjD graph II). 

The presence of such a regime leads to the thought of the existence of quasi-periodic 
regimes. To look for such a regime let us plot a bifurcation diagram of Poincare for devel¬ 
opment of basic limit cycle in case of increasing parameter a (FigJTb graph I). 

Another bifurcation, leading to the appearance of the toroidal surface, has been discovered 
in this system. An intersection of the toroidal attractor with the plane 1/3 = 0 forms a closed 
curve, shown in FigJH^. A further increase of parameter a causes the synchronization of tore 
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FIG. 7. a) Two-parametric bifurcation diagram in case of k = 2 (for other values of parameters 
and conventional symbols see FiglSl b) Poincare bifurcation diagram for development of the torus 
in case of 7 = 1, x = 10) C" = —2.8, r = 1.1, k = 0.9, h = 3.2, S = 0.4, ^ = 0.35, n = m = 3.2 
and increasing a, where graph I is the basic limit cycle, graph II - complicated periodic trajectory 
with separated region of attraction. 

frequencies, and finally an abrupt increase of vibrations amplitude, which shows the creation 
of a crucial new dynamical behavior. To clarify the character of the produced regime, let 
us analyze the Poincare section for the case of a = 14.6 (FigJHJo). The plotted cross-section 
is specihc for chaotic attractor, which provides reasons for statements on the existence of 
bifurcation of a quasi-periodic regime with a producing chaotic attractor. 

It turned out that system fl22|) provides another type of chaotic attractor creation, namely, 
intermittency. Let us fix 7 = 1, y = 50, G = —1.5, r = 0.1, k = 1.9, a = 0.2, h = 0.9, 

e = 0.18. 

Plotting the Poincare bifurcation diagram (Fig. |9^), we see that a limit cycle undergoes 
several period doubling bifurcations resulting in the chaotic attractor creation. But the 
development of chaotic attractor is interrupted suddenly and new complicated periodic tra¬ 
jectory appears which bifurcates in chaotic attractor as well at increasing n. Considering the 
hereditary sequences (Figj9b) for chaotic trajectories, we found that the graph of the map 
IFj+i = f{Wi) is close to the bissectrice at n = 4.25. As in the case with the Lorentz system, 
existence of narrow passage leads to the alternation of the chaotic and regular behavior of 
the system trajectories. 
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a b 

FIG. 8. a) The Poincare cross-section of the torus surface in case of a = 14 b) The Poincare 
cross-section of a chaotic attractor in case of cr = 14.6. Fixed parameters 7 = 1, x = 10) C' = —2.8, 
r = 1.1, K = 0.9, h = 3.2, <5 = 0.4, ^ = 0.35, n = m = 3.2. 



FIG. 9. a) The bifurcation diagram at increasing n. b) The graph of dependence ITj+i vs Wi 
at n = 4.25. The fixed values of parameters 7 = 1.49, x = 50, C* = —1.5, r = 0.1, k = 1.9, a = 
0.2, h = 0.9, C = 0.18, 5 = 0,8. 

VII. CONCLUSIONS 

Finally, we have studied the hierarchical sequences of the mathematical models for non- 
equilibrium media. Analyzing the wave helds in such media we have shown that derived 
models possess the wide set of localized wave regimes. In particular, the models with relax- 
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ation admit the periodic, multiperiodic, chaotic solutions. Spatially nonlocal models have 
in addition quasiperiodic and solitary wave solutions. All the models demonstrate the most 
bifurcations and scenarios of chaotic regimes creation. 

From the other hand, identifying internal variables with parameters undergoing fluctu¬ 
ations, one can consider these investigations as the problem on the dissipative strnctnres 
creation nnder the influence of noise. 
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